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Previously, the near-threshold pp — ► ppn reaction was studied with the use of transition operators 
derived from chiral perturbation theory (%PT) and the nuclear wave functions generated by high- 
precision phenomenological potentials. A conceptual problem in that approach was that the transition 
amplitude receives contributions from very high momentum components (above the cutoff scale of 
xPT) in the nuclear wave functions. In the present work, we avoid this problem by replacing the 
"bare" phenomenological potentials with Vi w-k> which is an effective potential derived from a bare 
potential by integrating out momentum components higher than a specified cutoff scale. The use of 
Viow-k is found to give an enhancement of the pp — ► ppir cross sections over the values obtained 
with bare potentials. Although this enhancement brings the calculated cross sections closer to the 
experimental values, the incident-energy dependence of the cross section is not well reproduced, a 
problem that seems to indicate the necessity of including higher chiral order terms than considered in 
the present work. 
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1 Introduction 



There have been many theoretical investigations PP- devoted to explaining the high-precision data 
for the total cross section of the near-threshold pp — > ppn° reaction |14l I15j . The initial surprise 
was that the measured cross section was larger than the values expected from the earlier calculations 
|16lll7j by a factor of ~ 5. Calculations in the phenomenological one-boson-exchange model indicated 
that heavy-meson (a and to) exchange contributions could account for the unexpectedly large cross 
section for pp — ► ppn The importance of heavy-meson exchanges in n° production is to be 
contrasted with their much less pronounced role in the charged-pion production process, which is 
dominated by the one-pion-exchange diagrams. Effective field theory (EFT), or more specifically, 
chiral perturbation theory (xPT) offers a systematic framework for describing the NN — * NNn 
processes at low energies. The leading-order term in xPT (the Weinberg- Tomozawa term) contributes 
to charged-pion production but not to n° production. xPT allows us to keep track of the contributions 
of higher chiral-order terms to the low-energy NN — ► NNn reactions jlJ|Sj. A point to be kept in mind, 
however, is that the NN — * NNir processes involve rather large momentum transfers, p ~ ^Jm^rn^ 
(ttt-tv = nucleon mass, = pion mass) even at threshold, and that this feature leads to the relatively 
slow convergence of the xPT expansion 6 j. The existence of the additional scale p ~ yjmfqm^ in 
the AW — ► NNtt reaction led Cohen et al.jHl to propose a new counting scheme in which the 
expansion parameter is \ = \Z m n/n^N instead of \w = rn^/miy employed in the usual Weinberg 
counting scheme. A thorough discussion on this and related topics as well as an extensive list of 
references can be found in a recent review by Hanhart ,13 . 

To maintain formal consistency in the xPT calculation of an inelastic nuclear process, one should 
derive from the same effective Lagrangian the relevant transition operators and the wave functions 
for the initial and final nuclear states. This type of calculation, however, has not yet been carried 
out. A practical and, in many cases, very useful method is a hybrid xPT approach |18|-|22|. in which 
the transition operators are derived from xPT but the nuclear wave functions are generated with 
the use of a modern high-precision phenomenological N-N potential. Hybrid xPT was applied to the 
pp — > ppn reaction in Refs. |U"!E1- These studies indicated: (1) There is a substantial cancellation 
between the one-body impulse approximation (IA) term and the two-body contributions, resulting in 
a cross section that is much smaller than the experimental value; (2) This feature seems reasonably 
stable against the different choices of phenomenological NN-potentials. 1 

A conceptual problem one encounters in these hybrid xPT calculations is that, whereas the transition 
operators are derived using xPT with the assumption that relevant momenta are sufficiently small 
compared with the chiral scale A x (p <C A x ~ 1 GeV) , the wave functions generated by a phenomeno- 
logical N-N potential can in principle contain momenta of any magnitude. A numerical calculation in 
Ref. 7 indicates that the transition amplitude receives non-negligible contributions from momentum 
components well above A x , a feature that jeopardizes the applicability of xPT. 

In a version of hybrid xPT called EFT* or MEEFT |19j-|21j. the contribution of the dangerously 
high momentum components in the wave functions are suppressed by attaching a momentum cutoff 
factor to the transition operators derived from xPT- 2 EFT* has proved to be extremely useful in 

1 Ref. 7 reports that the two different N-N potentials, Argonne V18 1231 and Reid soft-core |24| . give almost the same 
results for the pp — > ppn cross sections. 

2 Another important aspect of EFT* is that the low-energy constants appearing in the theory are constrained by the 
experimental data for the observables involving neighboring nuclei. This aspect of EFT*, however, will not be discussed 
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explaining and predicting many important observables for electroweak processes in few-nucleon sys- 
tems. Another possible way to suppress the contributions of high momentum components in hybrid 
%PT calculations is to attach a momentum cutoff factor to the wave functions. 3 Meanwhile, a sys- 
tematic method was developed by the Stony Brook group and others |251 126) to construct from a 
phenomenological NN-potential an effective NN potential that resides in a model space which only con- 
tains momentum components below a specified cutoff scale A. This effective potential, referred to as 
Viow-k; is obtained by integrating out momentum components higher than A from a phenomenological 
NN-potential, which in this context may be regarded as an underlying "bare" potential that resides 
in full momentum space. Vi ow -k represents a renormalization-group-improved effective interaction of 
a bare NN interaction. It has been found that, for a choice of A ~ 2 fm _1 , Vi ow _k reproduces low- 
energy observables such as the phase shifts (for p < A) and the deuteron binding energy with accuracy 
comparable to that achieved with the use of bare high-precision phenomenological potentials Fur- 
thermore, for any choice of bare NN-interactions (belonging to the category of modern high-precision 
phenomenological potentials), it has been found that the corresponding Vi ow -k generates practically 
the same half-off-shell T-matrix elements for p < A. This means that the low-momentum behavior 
(p < A) of the two-nucleon wave functions calculated from Vi ow -k is essentially model-independent. 

These developments motivate us to carry out a hybrid xPT calculation of the near-threshold pp — ► ppn 
reaction with the use of Vi ow _k. This type of calculation will substantially reduce the severity of the 
conceptual problem of momentum component mismatching that existed in the previous hybrid %PT 
calculations. It will thus allow us to examine more directly whether the transition operators derived 
from xPT up to a given chiral order are adequate or not. Furthermore, comparison of the results of a 
calculation based on Vi ow -k with those based on bare NN-interactions will also give information about 
the influences of the short-distance behavior of the NN-interactions on the NN — > NNtt reactions. In 
this context, it is informative to gather more examples of calculations that use bare NN-interactions. 
Therefore, in addition to a calculation based on Vi ow __k> we extend here our previous bare-potential- 
based calculations (carried out for the Argonne V18 and Reid soft-core potentials) to the Bonn-B 
potential [31] and the CD Bonn potential |33| . 

After describing the primary motivations of our work, we must mention that our present study is 
basically of exploratory nature and falls short of addressing a number of issues that warrant detailed 
studies. For one thing, we limit ourselves to the use of the Weinberg counting scheme, although 
it is important to examine the consequences of the counting scheme of Refs. As f° r the 

employment of Vi ow -k in hybrid xPT, there is a problem of formal consistency in that, whereas the 
relevant transition operators are derived in the framework of the dimensional regularization [IJ [7], 
Viow-k is based on the momentum cutoff scheme. Furthermore, the difference between the cutoff 
scales appearing in Ref. [7j and in V^-^ needs to be addressed. We relegate the study of these points 
to future work and concentrate on the examination of the consequences of the use of Vi ow _k in the 
present limited context; on the last point, however, we will give some brief remarks later in the text. 

The organization of this paper is as follows. Sect. 2 gives a brief recapitulation of the general framework 
of hybrid xPT, while we explain in Sect. 3 some technical aspects of numerical calculations we need to 
address as we work with Vi ow _k instead of the bare potential. The numerical results are presented in 

here. 

3 Insofar as the use of a momentum cutoff factor can be identified with the introduction of a projection operator onto 
a model space with a limited momentum range, applying the cutoff factor to the transition operators is equivalent to 
applying it to the wave functions. 
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Sect. 4 and compared with the data. Finally, Sect. 5 is dedicated to discussion and summary. 



2 Calculational Framework 

The formalism to be used here is basically the same as in Refs. [*fj [7j except that, for a calculation 
with Vi ow _k, some modifications (essentially of technical nature) are needed. Therefore, as far as 
the general framework of our approach is concerned, we only give a brief recapitulation, referring to 
Refs. HEj for details. 

2.1 Transition operators 

As in Refs. |IJ [7j, we derive the transition operators for the pp — ► ppir reaction using the heavy- 
fermion formalism (HFF) [*"7| of %PT based on the Weinberg counting rules. The relevant lagrangian 
is written as 

£ ch = £(°) + £« + £( 2 ) + --- . (1) 

Here C^ v > [p = 0, 1,2. . .) contains terms of chiral order v with v = d + in/2) — 2, where n is the 
number of fermion lines involved in the vertex and d is the number of derivatives or powers of m n . 
For our present study we only need the terms of v = and v = 1, which are given as follows: 

f 2 

£ (0) = i^Tr[d^d^U + m 2 n (U^ + U- 2)] + N{iv ■ D + g A S ■ u)N (2) 

£ (i) = -^-N{S-D,vu}N + 2cimlNNTr(U + U j -2) 
2totv 

+(c 2 --^-)N(v-u) 2 N + c 3 Nu-uN . (3) 

8m, jv 

Here U{x) is an SU(2) matrix that is non-linearly related to the pion field and that has standard 
chiral transformation properties; we use U(x) = \/l — [tt(x)/ f^] 2 + 2T-7r(x)// 7r . N(x) denotes the large 
component of the heavy-fermion field; the four-velocity parameter u„ is chosen to be = (1,0,0,0). 
Dfj,N is the covariant derivative, is the covariant spin operator, and = i[^d^ — '], where 
^ = y/U (x) The pion decay constant is taken to be f n = 93 MeV, and gA = 1.25. The values of 
the low-energy constants (LECs), c±,c<2 and C3, are given in, e.g., Refs. [4*| l2*%l [2*9*]: 

Cl = -0.87 ±0.11 GeV -1 , c 2 = 3.34 ± 0.27 GeV" 1 , c 3 = -5.25 ± 0.22 GeV" 1 . (4) 



The chiral order index v of a Feynman diagram is defined by v = 4 — Efq — 2C + 2L + ^ z/j, where 
En is the number of nucleons in the Feynman diagram, L the number of loops, C the number of 
disconnected parts in the diagram, and the sum runs over all the vertices in the Feynman graph |30j . 
We are using here the Weinberg counting scheme (HOj, although, as mentioned in the introduction, 
there exists a different counting scheme tailored to keep track of high-momentum flows involved in 
NN — * NNn reactions jHHEj. Furthermore, we limit ourselves to the tree- level diagrams and examine 
the consequences of employing Vi ow _k i n evaluating the transition matrix elements corresponding to 
these tree diagrams. 

The kinematic variables for the pp — > ppir reaction we use in this work and the relevant Feynman 
diagrams are shown in Fig. 1. As discussed in Refs. "HE], the lagrangian in Eq. (J""} leads to the 
transition operator 

J- — _|_ = 7-(Imp) _|_ ^-(Resc) ^ ^ 
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Figure 1: Impulse term (a) and rescattering term (b) for the pp — ► ppir reaction. In the text, the space 
components of the initial four-momenta pi and p2 in the center-of-mass system are denoted by pi and —pi, 
respectively; similarly, the space components of the final four-momenta p[ and p' 2 in the center-of-mass system 
are denoted by pf — q/2 and — pf — q/2, respectively. 



where represents the contribution of chiral order v. 7~( ^ comes from the one-body impulse 
approximation (IA) diagram [Fig. 1(a)] and is given by 



(2ir)3/2 ^5^2/,^' J 2m K 



7~( +1 ), which arises from the two-body rescattering diagram [Fig. 1(b)], is given by 



T (+l) = T (Resc) = _Z! V- , < fi- k 3 T 3 /-N 



where pj and p^' (j = 1, 2) denote the initial and final momenta of the j-th proton. The four- momentum 
of the exchanged pion is defined by the nucleon four-momenta at the ttNN vertex as kj = pj — pj, 
where pj = (E Pj , pj) , pj = (E p ?,pj) with E p = (p 2 + m 2 N ) 1 ^ 2 . The four-momenta of the final pion is 

q = (uj q , q), with uj q = (q 2 + m 2 /) 1 / 2 . The rescattering vertex function n{k, q) of Eq.(J7J) is calculated 
from Eq.©: 

^ [2ci - (c 2 - 9%A ) U;gfc ° - c 3 — ^ 
f 2 1 2 8mAr m£ 3 m 2 



«(*,,) = ^[2.*- (c 2 --^-)^-C3^] > (8) 



where /c = (fco, &) represents the four- momenta of the exchanged pion. 



2.2 Transition amplitude and nuclear wave functions 

We write the transition amplitude for the pp — > ppn reaction as 

T = {* f \T\*i), (9) 

where (!<&/)) is the initial (final) two- nucleon state distorted by the initial-state (final-state) inter- 
action. As briefly discussed in the introduction, in a formally consistent nuclear xPT calculation, the 
transition operator T and the N-N interactions that generate |<3?j) and \&f) are to be calculated to the 
same chiral order v from the common lagrangian. In hybrid xPT, we instead use a phenomeno- 
logical N-N potential to generate |<3?i) and In the present treatment this phenomenological N-N 

potential can be either Vi ow -k or a bare NN-interaction (see the introduction). 
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As described in Ref. |7j, we can apply the standard partial- wave decomposition to Eq. © and rewrite 
it into 



(X^, ms ,,m s ,?l r l^i si> m s2 ) E E ys f L? (py,™ s[ ,m s ^y s i L ;(pi,m sl ,m S2 ) 

1 2 " SfLfJfMfSiLiJiMi 



X 



E KmM (pflLfS^JfMflT^iq^piiLiSUMi) , (10) 
Here y^ff is the spin-angular function of the antisymmetrized two-proton state 

y J SL = 1 + { ~T +S E i L eM^LS)j}Yl ML (p)(^m Sl m S2 \SM s }(LSM L M s \JM}, (11) 

v M S M L 

where 5(ls)J 1S the NN scattering phase shift in the eigen-channel defined by the orbital angular 
momentum L, total spin S, and the total angular momentum J. l n denotes the angular momentum 
of the outgoing pion. It is convenient to introduce the reduced matrix element using the standard 
convention: 

(pflLfS^JfMflT^iq^piiLiS^JiMi) (12) 
- {-l) J f- M i ( _ J [ If £ * ) (Pf[L f S f ]J f || T Wm M II PilLiS^) . (13) 

Corresponding to the decomposition T = 7~ Imp + 7~ Resc in Eq.©, the reduced matrix element has 
two terms 

<Pf[L f Sf]Jf H^O?) \\pi[LiSi}Ji> 

=<p f [L f S f ]J f \\T^(q)\\p i [L i S i ]J i > + < Pf [L f S f ]J f \\T^ sc (q)\\ Pl [L l S i ]J l > • (14) 

Near threshold we can assume that the pp — > ppifi reaction is dominated by s-wave pion production, 
with the final pp states in the 1 Sq partial wave; this implies that we need only consider the 3 Pq partial 
wave for the initial pp state. With these constraints, the reduced matrix elements for the impulse and 
rescattering terms are given by 

^<Pf[%]\\T^ (q)\\M 3 Po]> = ~\ fff ^R 1So , Pf (p') 

4vr ,/(2vr) 3 2a; g /* J J 4?r ' 



x fi.( r #+-!?*-f')8(p , -?+?/2)R aFi)tPi ( P ), (15) 

771 JV 



' <Pf[ 1 So]\\T l % c (q)\\pd 3 Po]> = j=L [^Pr 1So , p (p>) 



K(k,q) p-k 
(2vr)3 k 2 -ml Po ' p * [P) ' 



(16) 



The radial functions, -R3p 0j Pi (p) and Ris , p f {p')i 111 Eas- lfTo)! and (fTB|) stand for the NN relative motion 
in the initial and final state, respectively. To obtain these radial functions, we first derive the JC-matrix 
by solving the Lippman-Schwinger equation in momentum space for a given NN potential; see Ref. |31j . 
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The calculated half-off-shell /T-matrix and N-N phase shift S(ls)J & ve the corresponding momentum 
space radial wave function as: 



R(LS)J,p (p) = i L cos{5 {LS)J ) 



$(p-Po) +v k (ls)j(p,Po,W) 



pi pl/m N -p 2 /m N 



(17) 



Here V means taking the principal-value part of the two-nucleon propagator, and po is the on-shell 
momentum defined by W = 2E Po . We note that the on-shell X-matrix is related to the phase shift as 
tan(5 {LS)J ) = -7rp m N K {LS) j(p ,p )/2. 

The choice of the four-momentum k of the exchanged pion has been a subject of investigations in the 
literature, see e.g., Refs 0J|32]. In a simple prescription, which came to be known as fixed kinematics 
approximation (FKA) 0j , we identify the four-momentums of the intermediate nucleon lines with the 
corresponding asymptotic values (ignoring thereby energy-momentum transfers due to the initial and 
final state interactions) and "freeze" all the kinematic variables at their threshold values. Thus FKA 
consists in using, in Eqs.Q and l)16[). k = (ko,k) = (m 7r /2,/j), k =p — p', and q = (u> q ,q) = (m^O). 
Meanwhile, in Ref. [7j, k was chosen in such a manner that energy- momentum conservation at each 
vertex in Fig.l(b) should be satisfied, i.e., 

k = {k , k) = (Ep - E f ,_ ?/2 ,p- p' + q/2) , (18) 

where p (p') is the relative momentum of the two protons before (after) pion emission, defined in 
Fig. das: 

Pi = -P2=P, Pi=P Pi=-P -g- v 19 ) 

A schematic study in Ref. |32j indicates that, when the final-state NN interaction is included in the 
rescattering diagram, FKA is an appropriate choice, but that, when initial-state interaction is included, 
the situation is more complex. Thus the choice of k is still an open issue. In the present work therefore 
we give numerical results for the choice of k given in Eqs. (|18j) and (|19|). as well as for FKA. The bulk 
of our calculation will be done with the use of k given in Eas.(|18|) and ([190: the results corresponding 
to FKA will be presented with due remarks attached to them. 



2.3 Cross sections 

The total cross section at energy W(= 2Ep i ) in the center-of-mass frame for the reaction pp — ► ppir 
is given by [7j 

Clir) 4 E rim , 

LiSiJiLfSf J flir 

(20) 

where E Pf = {(W — oj q ) 2 — g 2 } 1 / 2 /2, pf = \J E 2 f — m 2 N , and the maximum momentum, q m , of the pion 
is given by q m = \J {(W — 2rriN) 2 — m1.}{(W + 2mN) 2 — m 2 } / 4W 2 . Here pi (= \pi\) is the asymptotic 
relative momentum of the initial pp states and E^ t = ^Jp 2 + m? N . Since we have already specialized 
ourselves in the threshold pion production, we need not deal with the general expression in Ea.(|20|): 
we can limit [LfSf]Jf to and [LiSi]J{ to 3 Pq- 
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2.4 NN interactions 



As discussed in the introduction, the main purpose of the present work is to carry out a hybrid 
calculation of the pp — ► ppn reaction with the use of Vi ow _k, which resides in reduced Hilbert space 
characterized by the constraint p < A. Specifically, we use here the Vi ow _k derived from the CD-Bonn 
potential [SI- 4 We are also interested in comparing the the results of this calculation with those 
of hybrid xPT calculations based on standard high-precision phenomenological potentials, which we 
refer to as "bare" interactions. Regarding a bare NN potential case, in order to augment the examples 
given in Ref. 0, we shall carry out additional calculations with the use of the Bonn-B potential and 
the CD Bonn potential. 



3 Numerical Calculation 

The numerical evaluation of the scattering amplitude and the cross section follows closely the method 
employed in Ref. [7]. The calculation of the rescattering amplitude [Eq. (|16|)] can be readily done 
in momentum (p)-space, whereas it is technically easier to carry out a numerical evaluation of the 
impulse amplitude [Eq.pS])] in coordinate (r)-space [7j. Since the calculational method for the case of 
a bare NN-potential was explained in detail in Ref. 0, we only describe here modifications that need 
to be made when we use Vi ow _k instead of a bare potential. 



3.1 Principal- value integral 

The principal value integral appearing in Eq. (|17j) is usually rendered amenable to numerical calculation 
in the following manner. If we need to numerically evaluate the integral 

lsP T dk 7^- (21) 

we may convert this expression into an ordinary integral by subtracting zero: 



q 2 — k 2 Jo q 2 — k 2 

OO 1 

In a calculation that involves Vi w-kj the upper limit of k-integration is a finite value (A), so that we 
encounter an integral like 

l A ^vj\k-l^, (23) 

where, for the sake of definiteness, we may assume q < A. Since V J Q A dk q (^li ^ 0, the procedure 
used for I needs to be modified as 

■ ' /(*) r A „ /(«) , ^ f A i, /(?) 



q 2 — k 2 

A 1 fid) A + q 



We are grateful to T.T.S. Kuo for allowing us to use a computer code to generate Vi ow _k developed by his group. 
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3.2 Calculation of the impulse-term amplitude in r-space. 

As mentioned, the numerical evaluation of the impulse-term amplitude can be done more conveniently 
in r-space than in p-space. In a case involving a bare potential, switching from the p-representation 
to the r-representation can be readily performed using the standard Bessel transformation, 



oo 

R{LS)J, P0 (r) = \j^i L P 2 dpjL(pr)R( LS) j, P0 (p) , (25) 
and the well-known identity, J °° j L(p r )j L(p r ')p 2 dp = -^z5(r — r'). The result is given by [7] 
-i <Pf[ 1 S ]\\T^ Q ( q )\\p l [ i P ]>= 1 f fdrr 2 R^ pf {r) 



(1 + ^-)Qh(qr/2) - ^Lj ( qr /2)(± + -)] Rsp^r) . (26) 
IrrtN mjv dr r J r 

The usefulness of this method, however, diminishes in the case of Vi ow _k, where the momentum 
integral does not run to oo but stops at A, and hence we cannot use the above-quoted orthogonality 
of the spherical Bessel functions: J Q 3L{pf)jL{pi' l )p 2 dp ^ ^zd( r ~ r ')- We therefore use the following 
procedure. In evaluating the impulse amplitude [Ea. (J15|) ] for Vi ow _k; we first integrate out the 5- 
function, and then divide the range of p-integration in two intervals as follows. 

J -<Pf[ 1 So]\\T^ ( q )\\ Pl [ 3 Po}>- ~ — 



x{/ dpp 2 dxRx S Ap) - ■ (-q+ — SL p) R 3p (I) 
Jo J-i u ' y} L< tun J 

+ [ A dpp 2 [ +1 dxRi Sos (p) \ l ~ • (-q+ ^Lp )} i?3 Poi n (I) } . (27) 

Jrt r J — 1 L t m N 



Here p = \p \, and x denotes the cosine of the angle between p and q, i.e., p '■ q = pqx; I = p + q/2, 
and the momentum p c is chosen to lie between pf and p m , where p m is the solution of the equation 
Pi = \Jpfn + PmQX + (? 2 /4 for a given value of x. The merit of dividing the p- integration range in the 
two intervals is that each p-space integral in Eg. 1)271) contains only one principal- value part coming 
from either the initial or the final NN relative-motion propagator. For instance, in the second integral 
of Eq. (|27j) the final state radial wave function takes the following simple form 

Ri So , Pf (p) = cos[5(}S )}m N Kls 2 oip ' P 2 f) , 

P f P 

where Pf(< p c ) is the final-state on-shell momentum. We evaluate the first term in Eq. (|27|) directly in 
p-space. The second integral in Ea. (|2*7|) is calculated using a modified Bessel transformation outlined 
in the appendix. 

Since the above-described method for carrying out the r-space calculation of the 1-body amplitude 
with Viow-k is somewhat involved, it seems safer to check its validity using some pilot calculation. 
If in deriving T we assume q = (the "<f = approximation"), the evaluation of the transition 
matrix element is drastically simplified, and even with Vi ow -k we need not resort to the above lengthy 
prescription. We therefore consider it informative to compare the results of calculations with and 
without the "<f = approximation". This comparison is given in Appendix B. 
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Figure 2: <r lmp , the pp — > pp7r° cross section calculated with the impulse term alone. The cross section is given 
in units of rj 2 , where the "velocity" 77 is defined as the maximum pion momentum at a given laboratory energy 
Ehab of the incident proton divided by the pion mass, ra^ = 135 MeV. 

4 Numerical Results 

We now present the pp — > ppir® cross sections calculated with the use of Vi ow _k [2S] as well as typical 
bare NN interactions. We consider here three examples of the bare potential: Argonne V18 |23| . 
Bonn-B |34j . and CD-Bonn [SHI- The results for Argonne V18 are taken from Ref. 0, while those for 
Bonn-B and CD-Bonn have been newly calculated in the present work. It is informative to study the 
individual behavior of the impulse- and rescattering-term contributions before discussing the behavior 
of their combined contributions. 

4.1 Contribution of the impulse-approximation term 

We first consider <r Imp , the total pp — > ppir cross section calculated with only the impulse-approximation 
amplitude retained; viz., in evaluating Ea, ()20j) . we replace <pj[LfSf] Jf \\ Ti^{q) \\pi\LiSj\ Ji > with 
<p/[ 1 S'o] || T\ =o(?) llPiP^o] >• Fig- El shows <j Imp as a function of i?Lab> the incident proton energy 
in the laboratory system. We note that, for the three representative bare NN-potentials, cr Imp varies 
up to 40 %. These variations are a measure of ambiguity inherent in a calculation that uses a bare 
NN potential. The fact that the short-distance behavior of bare NN potentials is not controlled with 
sufficient accuracy underlies this instability. Fig. El indicates that the use of Vi ow -k leads to a value 
of <r Imp that is significantly smaller (by a factor of 3 or more) than those for the bare potentials. A 
plausible explanation of this difference is as follows. In the one-body transition diagram [Fig.^a)], the 
large momentum transfer {p ~ ^/m„-mjv) between the two nucleons needs to be mediated by the N-N 
potential. 5 Now, by construction, Vi ow _k only contains momentum components below A = 2fm~ 1 , 
whereas the bare potentials carry very high momentum components (albeit in a rather arbitrary man- 
ner). We can expect that the absence of those high-momentum components in Vi ow -k suppresses the 
contribution of the one-body transition diagram. 

5 As pointed out in Ref.|7|, with the use of a bare NN potential, one needs to take the upper limit of p-integration very 
high (up to p ~ 2 GeV/c) before the integral starts to show a sign of convergence for both the impulse and rescattering 
amplitudes. 
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Figure 3: ,j Rcsc , the pp — > ppn cross sections calculated with the rescattering term alone. For the explanation 
of the quantity 77, see the caption of Fig.l. 

As mentioned, we are using in the present work the Vi ow __t derived from the CD Bonn potential. It is 
known, however, that, so long as one starts from a bare potential that belongs to the category of mod- 
ern high-precision phenomenological potentials, the resulting Vi ow _k is practically model-independent 
in the sense that the half-on-shell .fT-matrices corresponding to different bare potentials are nearly 
indistinguishable [251 I26| . This means that <r Imp calculated with Vi ow _k corresponding to any realistic 
bare potential would lie close to the solid line in Fig. |2j Thus the use of Vi ow -k results in a significant 
reduction of model dependence in our calculation. 

4.2 Contribution of the rescattering term 

Fig. |21 gives fj Resc , the total pp — > pptr cross section calculated with only the rescatttering term 
contribution retained; i.e., in evaluating Eq. (|2fl|) . we replace <pf[LfSf]Jf \\Ti^(q) \ \pi[LiSj\Ji> with 
<pj[ 1 So] ||^ R Jfo(<z) ||Pi[ 3 ^o]>- The figure indicates that, for the three different choices of the bare 
NN-potential, cj Resc shows variations of about 30 %, while the use of Vi ow -k leads to a Resc that lies 
more or less within the range of these variations. In the rescattering diagram [Fig. 1(b)], a substantial 
fraction of the momentum transfer between the two nucleons can be carried by the exchanged pion, 
and therefore the NN interactions need not directly support a large momentum transfer. This feature 
explains why the change in <j Resc is less pronounced than a Imp as we switch from the bare potentials 
to Vi ow _k. It is worth re-emphasizing here that, although <7 Resc (Vk> w _k) i n Fig. |3] was obtained with 
the Vi ow _k derived from the CD-Bonn potential, the result should be considered model-independent 
in the sense discussed in the preceding subsection. 

4.3 Combined contributions of the impulse-approximation and rescattering terms 

We now consider the total pp — ► ppw cross section, a, calculated with the full transition amplitude 
consisting of the one-body and two-body terms; thus a is obtained from Eq. (|2()(l with the transition 
amplitude given by Eqs. (|14|). (|15|) . (|16|) . The results are shown in Fig. 0] for the three choices of 
the bare potential and for Vi ow -k> along with the experimental values of cr pp ^ pp7T o. Fig. 0] indicates 
that the use of Vi ow _k leads to a rather visible enhancement of a over the results obtained with the 
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Figure 4: The total cross sections calculated from the full amplitude consisting of the impulse and rescattering 
terms. The experimental data points are taken from Ref. \TQ (solid circles) and Ref. ^5] (solid squares). For 
the explanation of the quantity rj, see the caption of Fig.l. 

bare potentials. This enhancement is related to the suppression of the impulse-approximation ampli- 
tude corresponding to Vi ow -k- As pointed out in the earlier calculations jUOIl, the impulse and 
rescattering amplitudes tend to interfere destructively, and in the case of a bare NN interaction the can- 
cellation between the two amplitudes is quite substantial, leading to a significantly suppressed value of 
a as compared with the individual magnitudes of <r Imp and cr Resc . The smaller impulse-approximation 
amplitude obtained with the use of Vi ow -k somewhat diminishes the extent of this destructive inter- 
ference, resulting in a larger value of a. The enhancement of the cross section obtained with Vi ow -k 
brings the calculated values of a closer to the experimental values. It is to be noted, however, that the 
energy dependence of a obtained with Vi ow _k differs significantly from the experimentally observed 
behavior. We remark once again that the o"(Vi ow -k) shown in Fig. 0]should be essentially independent 
of the choice of a bare potential from which Vi ow _k is derived (see subsection 4.1). 

The above results correspond to the case where the four-momentum k of the exchanged pion is specified 
according to Eqs. (|18j) and (|19|) , As discussed earlier, however, there is an argument that favors FKA 
in a certain context [32] ■ It therefore seems informative to repeat our calculation with the use of FKA. 
The a corresponding to this case is shown in Fig. 03 We observe that the cross sections obtained in 
FKA are smaller than those in Fig. |IJ representing a larger deviation from the experimental values, 
and that the energy dependence of a does not resemble the experimentally observed behavior. 

However, we need to discuss here the dependence of our results on the values of the LEC, c\, C2 and 
C3. The above results were obtained for the "standard values" of ci, C2 and C3 given in Eq.jjlJ. These 
were originally deduced in Ref. |2B] and quoted in Ref. [7j as "parameter set I" . The allowed ranges 
of these LECs were discussed in Ref. 0122]) where, in addition to the parameter set I, two more sets 
were considered as examples of other possible choices. For convenience, we tabulate these three sets 
of parameters: 

Parameter set I 
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Figure 5: The total cross sections calculated from the full amplitude consisting of the impulse and rescattering 
terms evaluated in the fixed kinematics approximation. The experimental data points are taken from Ref. |14| 
(solid circles) and Ref. Jl] (solid squares). For the explanation of the quantity 77, see the caption of Fig.l. 

Cl = -0.87 GeV -1 , c 2 = 3.34 GeV" 1 , c 3 = -5.25 GeV" 1 (28) 
Parameter set II 

Cl = -0.87 GeV" 1 , c 2 = 4.5GeV-\ c 3 = -5.25 GeV" 1 (29) 
Parameter set III 

cx = -0.98 GeV -1 , c 2 = 3.34 GeV -1 , c 3 = -5.25 GeV -1 . (30) 

To get a measure of the sensitivity to the choice of the LECs, we have repeated our calculation of a 
for Vi ow „k using the parameter sets II and III. The results are shown in Fig. El together with those for 
the set I; in fact, since the set II gives practically the same result as the set I, we give in the figure 
only the results for the sets I and III. Fig. El indicates that the set III, which differs from the set I 
only by a modest 12% change in c\ = —0.98 GeV^ 1 , enhances a considerably, bringing the calculated 
values of a closer to the experimental values. However, the energy dependence of the theoretical a 
remains dissimilar to the experimentally observed behavior. Figs. 01 El and El seem to suggest that, in 
order to fully explain the magnitude and incident-energy dependence of the pp — > ppw cross sections 
near threshold, one probably needs to include terms of chiral orders higher than considered here. 
We remark in this connection that the possible importance of two-pion exchange diagrams in a xPT 
calculation of the pp — > ppir reaction was pointed out in Refs. |35[ I36j . 



5 Discussion and Summary 

We have carried out a hybrid xPT calculation of the cross section a for the s-wave pion production 
reaction, pp — ► ppn , with the use of Vi ow _k- Vi ow _k is a low-energy effective potential derived from 
a high-precision phenomenological potential (called a "bare" potential in our context) by integrating 
out momentum components higher than A ~ 2 fm _1 . The results obtained with Viow-k are compared 
with those obtained with the three representatives bare potentials, ANL V18, Bonn-B and CD-Bonn. 
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Figure 6: The total cross section a calculated with the full amplitude (including both the one-body and 
two-body amplitudes) for Vi ow _k- The results obtained with the parameter sets I and III are compared. The 
experimental data points are taken from Ref. (solid circle) and Ref. ^3 (solid square). 

The principal features of our calculation based on Vi ow -k are summarized as follows. 

(1) A hybrid xPT calculation based on a bare potential has the "momentum mismatch" problem that 
the initial and final nuclear wave functions generated by the bare potential involve very high momenta, 
whereas the transition operators derived from %PT can be used only within a limited momentum range 
(p < 300 MeV). This momentum mismatch problem is significantly mitigated by the use of Vi ow -k- 

(2) Reflecting the fact that the short-distance behavior of the bare potential is not well controlled, 
the cr's calculated with the above-mentioned three bare potentials exhibit ~40 % variance. This kind 
of model dependence practically disappears with the use of Vi ow _k, since different choices of a bare 
potential lead to practically equivalent Vi ow -k's [2E]- This feature allows us to better focus on the 
question whether the transition operator for the pp — > ppir reaction derived from xPT up to next-to- 
leading order is adequate or not. 

(3) The calculation with Vi ow -k enhances a over the values obtained with the bare potentials, and, 
with certain choices of the relevant LECs, a can come close to the experimental values for some range 
of the incident energy. It is however unlikely that the magnitude and energy dependence of a can be 
fully reproduced with the transition operators considered in this work. Thus it seems necessary to 
consider higher-order transition operators. 

For formal consistency, it is desirable to go beyond the hybrid approach by employing N-N 

potentials derived from xPT. This is however a major task relegated to the future. We remark in this 
connection that, for reactions that only involve the rearrangement of the nucleons, there has been much 
progress in constructing a framework that is formally consistent with effective field theory |37l I38j . 

A related issue is that we concentrated here on the consequences of changing the nuclear wave functions 
from those generated by the bare NN potentials to those generated by Vi ow _k, without taking account 
of the possible renormalization of the transition operators due to the truncation of model space. 6 As 
is well known, a reduction of nuclear model space in general entails a corresponding modification of 

For discussion of some formal aspects of the use of Vi ow -k in hybrid x?T, see Ref. PSI- 
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operators for the nuclear observables. Again, to fully address this issue, we need to go beyond the 
hybrid used in this work. In the present study of exploratory nature, we have concentrated on 
a hybrid xPT evaluation of the transition operators that arise from the tree diagrams. At this chiral 
order there are no loop corrections to the LECs (cj, i = 1, 2, 3) and the other coefficients appearing 
in C c h like, e.g. g^. We may therefore expect that, although the scale of xPT (A x ~ 1 GeV) is larger 
than the momentum cutoff scale (A ~ 2 fm _1 ) used in deriving Vi ow _k, this difference does not lead 
to a drastic renormalization of the transition operators. To turn around the argument, the issue of 
the renormalization of the transition operators is coupled with the treatment of higher chiral-order 
terms, and these two aspects need to be addressed simultaneously. This important question, however, 
is beyond the scope of our present study, which is limited to the tree diagrams. 

As an immediate follow-up of the work described here, we are studying j4()j the expected important 
contributions from the two-pion exchange diagrams [351 156 in a hybrid xPT calculation with Vi ow _k. 
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Appendix A 

With the use of Vi ow _k, the numerical evaluation of the second integral of the impulse amplitude 
[Eq. (j27jl ] requires the "modified" Bessel transformation, 



R(LS)J, P0 ( r ) = ^J~ iL J P 2d PjL(P r )R(LS)J, P0 (p) ■ (31) 

However, since jL(p r )jL(p r ')p 2 dp ^ ^pzd( r ~ r ')> the inverse transformation for Eq. ()31|) gets com- 
plicated, which presents us from arriving at an r-space expression similar to Ea, ()26|) . We therefore 
replace R(ls)J,p (p) i n Eq. (|3*Tj) with the K-matrix expression Eq. lfTTJ) and then we use in Eq. (|27j) the 
Bessel transformation 



oo 



R(LS)J, Po (k) = y-» / r drj L (kr)R iLS)J:P0 (r), (32) 

where k < A is understood. In our numerical evaluation of the impulse amplitude with Vi ow -k) only 
the initial wave function in the second integral of Ea. ()27|) is evaluated using the above prescription. 



Appendix B 

To check the validity of the numerical techniques described in the text, it is informative to compare the 
<7 Imp calculated for Vi ow -k using the prescription explained in subsection 3.2 with the cr Imp obtained 
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Figure 7: <r Imp calculated for Vi ow -k without or with the q — approximation. The solid (dotted) line 
corresponds to the calculation without (with) q = approximation. For the explanation of the quantity 77, see 
the caption of Fig.l. 

in the u q = approximation", wherein the transition operators are derived under the simplifying 
assumption that the out-going pion has no momentum, q = 0. [In evaluating the phase space, we do 
treat q as a variable.] In the q = approximation, the impulse term can be evaluated in p-space in a 
straightforward manner without any complicated handling of the principal value integrals. The results 
for the two cases, with and without the q = approximation, are shown in Fig. [7J and we compare this 
figure with Fig. 2 in Ref. [7J, which presents <7 Imp calculated for the bare potentials with and without 
the q = approximation. According to Fig. 2 in Ref. [7J, <7 Imp obtained in the q = approximation 
tends to become somewhat larger than cx Imp obtained without using the q = approximation. The 
fact that Fig. [7J exhibits a similar tendency may be taken as an indication that the somewhat lengthy 
procedure we use in handling the principal-value integrals for the Vi ow -k case is reliable. 
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